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The end-to-end energy - energy correlations of random transverse-field quantum Ising spin chains 
are computed using a generalization of an asymptotically exact real-space renormalization group 
introduced previously. Away from the critical point, the average energy - energy correlations decay 
exponentially with a correlation length that is the same as that of the spin - spin correlations. 
The typical correlations, however, decay exponentially with a characteristic length proportional 
to the square root of the primary correlation length. At the quantum critical point, the average 

-const-L I ^ wnereas the typical correlations decay 
e~ K ' , with K a random variable with a universal distribution. The critical energy- 



correlations decay sub-exponentially as Cl 
faster, as 

energy correlations behave very similarly to the smallest gap, computed previously; this is explained 
in terms of the RG flow and the excitation structure of the chain. In order to obtain the energy 
correlations, an extension of the previously used methods was needed; here this was carried out via 
RG transformations that involve a sequence of unitary transformations. 



I. INTRODUCTION 

The random transverse field Ising model is the proto- 
typical model of a quantum system with quenched ran- 
domness. It undergoes a novel quantum phase transition 
at zero-temperature, and over a wide range of parame- 
ters, exhibits unusual low temperature behavior. In one 
dimension, many of the low energy properties have been 
found exactly, initially by McCoy and Wu 1,2 and in much 
detail recently. 3-5 The behavior in higher dimensions has 
been shown to be in the same general class as the one 
dimensional system. 6 ' 7 

In addition to its theoretical interest, models with 
similar behavior have been argued to be relevant for 
the low temperature properties of heavy fermion materi- 
als, with randomness and proximity to a quantum criti- 
cal point playing key roles in producing non-fermi-liquid 
behavior. 8 

In this paper we study the end-to-end energy and 
transverse spin correlations of long but finite random 
transverse field Ising chains. The average energy and 
transverse spin correlations are found to have rather dif- 
ferent behavior than the order parameter correlations. 
Computing them introduces new difficulties which com- 
pel us to further develop the general renormalization 
group (RG) structure. The formulation presented here 
should be useful for other random quantum systems. 

The organization of this paper is as follows; in the re- 
mainder of this section we review some aspects of the 
random Ising model and introduce the quantities of in- 
terest. In Sec. II we develop and apply a unitary- 
transformation approach to the real-space RG. In Sec. 
Ill and IV Laplace-transforms of the energy and trans- 
verse spin correlations are derived and used in Sec. V 
to obtain the average and typical correlations, and more 
general information about the distributions. Finally, Sec. 
VI presents conclusions. Some technical details are con- 
fined to an Appendix. 
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FIG. 1: The Hamiltonian of the transverse field Ising model. 
Each site is a spin- 1/2 that interacts via Ising exchange 
with its nearest neighbors and can be nipped by the local 
x-magnetic field. 



A. Random Transverse Field Ising Model 

The Hamiltonian of the random quantum Ising model 



is 



ft = - ( J «+i°f °f+i + hiaf) 



(1) 



with each site having two states, a z — ±1, with quantum 
fluctuations between them caused by the transverse, <7 X , 
fields. The system is illustrated in Fig. 1. Note that there 
are no magnetic fields in the z-direction, so that there is 
a global symmetry of inversion about the xy spin plane. 
The presence of z fields would break this symmetry and 
change the low energy physics radically. 

The quantum Ising model exhibits a quantum phase 
transition in its ground state when the nearest neigh- 
bor interaction and the transverse field are of compara- 
ble strength. In a non-random model this occurs when 
J = h. In a random system, where the J's and the h's are 
drawn independently from some distributions, the tran- 
sition occurs when log h = log J, where the over-bars 
denote averaging over the randomness. A convenient 
parametrization of the proximity to the transition is 



= log hi - log J i 
var(log h) + var(log J) 

with S positive in the disordered phase. 



(2) 
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B. Real Space RG 

A powerful route to analytic information on this sys- 
tem is a real space — or energy space — RG method 
that is a generalization developed by one of us, 4 of an 
RG introduced by Ma, Dasgupta and Hu. 9,10 The real 
space RG is carried out by decimating the term in the 
Hamiltonian — a site (ha x ), or a bond (Ja z a z ) — with 
the strongest interaction; second order perturbation the- 
ory results in new effective couplings. In the case of 
decimating a bond, a cluster forms with a renormalizcd 
transverse field; in the case of decimating a spin, a new 
Ising interaction that couples its two neighbors forms. 
The Hamiltonian preserves its form, with effective bonds 
coupling spin clusters, and the energy scale — the maxi- 
mum remaining coupling — reduced. The effective bond 
strengths and lengths, as well as the effective transverse 
fields on the clusters of spins and their moments can be 
computed. 

As the energy scale is systematically reduced, the dis- 
tributions of the effective couplings become very broad 
for small d. Concomitantly, the averages of many quan- 
tities in the ground state are determined by rare tails of 
their distribution. In particular, the average order pa- 
rameter correlations, 



C zz (n, n + r) 



< o~ z a. 



>, 



(3) 



at large separations r, arc dominated by the rare pairs 
of spins that are not decimated until they join together 
into the same spin cluster. 4 



C. Logarithmic Energy Scaling 

Many properties of the random quantum Ising model 
can be understood in terms of the scaling behavior of the 
cluster sizes, bond lengths, and coupling strengths with 
the energy scale and the deviation from criticality, 5. 

At the critical point, the distributions of hi and Ja+i 
become infinitely broad as the energy scale, fi, ap- 
proaches zero; the random quantum critical point is thus 
an infinite randomness fixed point. 6 ' 11,12 At this fixed 
point the distribution of cluster and bond lengths, the 
logarithms of the interactions in units of O, 



C^log^, 
/3 = log£, 



(4) 



and the deviation from criticality, S, all scale with the 
logarithm of the energy scale, 



r = io g 



0/ 



(5) 



Here fij is the initial energy scale set by the strongest 
couplings, and f2 is the magnitude of the largest remain- 
ing couplings after the stronger ones have been deci- 
mated. The scaling can equivalently be given in terms of 



a length scale, I — for example the length of an effective 
bond — the scaling of log-energies at fixed I is of the 
form 



C = zVl, 

r = 7 A 



(6) 



where 7 and z are scale invariant random variables. 

Various basic results follow directly from this scaling. 
In particular, the linear number density of remaining 
spin-clusters at scale T is 



no 

p2 ' 



(7) 



with no a non-universal prefactor inversely proportional 
to the original bond lengths. 

An example of the scaling of log-energies with length is 
the gap, Ei — Eq, between the ground state and the first 
excited state of long finite chains: Analytic and numeri- 
cal results show that near the critical point the logarithm 
of the gap is broadly distributed on the scale \[L. Indeed, 
the distribution of — \a(E\ — E )/^/L attains a univer- 
sal scaling form in the large L, small S limit. 3 From the 
RG structure, this can be seen by noting that the gap 
is the energy scale of the chain when it has only one re- 
maining cluster — and thus only one unfrozen degree of 
freedom. Therefore the gap is approximately ~ e~ ri 
with oc \[L the sample-specific scale at which this last 
cluster is decimated. 

In long finite chains of length L, the end-to-end spin 
correlations < o~q<j z l > are a useful probe of the long 
length scale ordering tendencies. The distributions and 
moments of these can be calculated exactly in the asymp- 
totic limit of long chains and small 5; these compare 
well with numerical results. 3 The distributions can be 
expressed usefully in terms of 



A 2 



■lGgC**(0,L) . 



(8) 



This logarithm of the correlations scales with \[L at the 
critical point, with a broad distribution on the same 
scale. The average correlations, however, decay much 
more slowly: only as C zz oc 



L ' 



D. Ordered and Disordered Phases 

When S is non-zero but small, there are two scaling 
regimes. At early stages of the decimation process, clus- 
ters and interactions are not "aware" of being non crit- 
ical. In this regime the critical scaling holds. At longer 
scales, however, there is a crossover to an off-critical 
regime. The crossover occurs when the typical cluster 
sizes and bond lengths are of order of the correlation 
length 



1 



(9) 



3 



and the log-energy scale is of order 

r> 



1 

6 ' 



(10) 



At scales larger than £, the behavior is characteristic of 
one of the two zero-temperature phases and thus depends 
on the sign of S. 

At low energies in the ordered and disordered phases, 
the scaling between energy and length scales is different 
from that at the critical point. For small 5, in both phases 



0~f 



-z(<5) 



with the effective dynamical exponent, 

1 



(11) 



(12) 



near the critical point. 

The distributions of the log-interactions also change 
form. In the disordered phase, the distribution of the 
effective fields does not continue to broaden and (3 ~ 
r x ~ j. But the effective bonds become longer and 
longer and weaker and weaker with the distribution of the 
In J's broadening rapidly, with typical £ ~ T. In the dis- 
ordered phase, the average order parameter correlations 
decay exponentially with the correlation length £ w ^. 
Nevertheless, the typical correlations decay much faster; 
for example, end-to-end correlations of almost all samples 
decay as e~ 2SL . 13 More precisely, as L — > oo, the distri- 
bution of of the scaled log-correlation function, A z /L, 
approaches a delta function peaked at A z /L = 25. The 
average correlations are thus dominated by exponentially 
rare samples that happen to have anomalously strong ex- 
changes and/or anomalously weak random fields. 

In the ordered phase, 6 < 0, at low energy scales the 
clusters become bigger and bigger, eventually encompass- 
ing the whole system. The transverse fields on these 
clusters — the gap between the symmetric and anti- 
symmetric combinations of their two "ordered" states 
- concomitantly continues to become more and more 
broadly distributed. The remaining bonds, on the other 
hand, stay relatively short and their distribution does not 
continue to broaden. These fields and bonds thus play 
opposite roles in the two phases; as discussed below, this 
is a general consequence of duality. 



E. Duality 

As for non-random classical and quantum Ising models, 
there is a dual description of the random quantum chain 
in terms of bond variables. Instead of using the states 
| |> 5 I l>> on eacn site , one can use the states of the 
bonds. This is done by assigning |+ >, |— > to the bond 
if the two spins surrounding the bonds are | |T>i I II > 
or I IT> 5 I II > respectively. These are domain- wall vari- 
ables. In the new Hilbert space, if we choose the quan- 
tization axis to be x rather than z, the Hamiltonian has 
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FIG. 2: A chain that terminates with the site is dual to 
a chain that terminates with a site 0',which experiences no 
transverse field. The bond operator —hoa^a^, is then the 
edge energy operator of the dual chain. 



the same form, but with h and J exchanged. The duality 
is summarized in the following table: 



S 

Jn n+1 

K 



-s 
K 

Jn n+1 



(13) 



o-„<t: 



n u n+l 



If the S dependence of the distributions of the random 
couplings has the form pj(J = X, 5) = pu{h = X, —6), 
the random model is self dual with 6 => —S. More gen- 
erally, it will not be exactly self dual. But from the 
definition of d in terms of pj and ph (Eq. 2), and the 
universality at low energy scales (in particular, of the 
distributions of the effective couplings), we expect that 
the asymptotic behavior at low energies and small S will 
indeed be self-dual for any well-behaved distributions. 

As we are interested in end-to-end correlations of finite 
chains, we must consider what happens to the ends of a 
chain under the duality transformation. Let — /ioO"o be 
the energy operator on the left end site. Under duality, 
this site will be mapped into a bond with corresponding 
energy operator —hoa^af, (see Fig. 2). The bond in the 
dual chain implies a new site 0' that corresponds to the 
domain wall variable to the left of the first spin in the 
original chain. This new site, 0', thus carries the informa- 
tion about an arbitrarily fixed boundary condition, e.g., 
| |>, with respect to which the leftmost domain wall, 
and hence the original end spin ctq, is defined; it thus 
cannot be entirely forgotten. But in the dual Hamilto- 
nian, the operator Oq, does not appear, so that the dual 
end transverse field is zero, and Ctq, is time independent. 
The same is true at the other end, where the extra spin 
is needed for the original cr£ to be defined. Superficially, 
the dual chain appears to have one extra degree of free- 
dom associated with each end. But the orientation of the 



4 



dual edge-spins is entirely a convention, and therefore the 
additional degrees of freedom have no effect. 

Note that in the special case in which the original chain 
starts with a spin to which no transverse field is applied, 
the duality transformation yields a chain with the first 
bond turning into a site. More generally, any site in the 
original chain on which there is no transverse field corre- 
sponds under duality to a break in the chain that divides 
it into two uncoupled parts. In the original variables, 
there are concomitantly two disconnected subspaces in 
which the spin that cannot flip has a z = ±1 respectively. 
The parts of the chain to the left and to the right of this 
spin are thus independent of eachother. 



F. Energy-Energy Correlations 

Our goal in this paper is to understand the energy den- 
sity (E-E) correlations of the random chain. These would 
be particularly interesting at non-equal times, as they 
would then yield information on the transport of energy 
which is the only locally conserved quantity in this sys- 
tem. Unfortunately, in the bulk of the chain correlations 
are very hard to calculate for reasons discussed in rcf. 
4. Therefore we study the somewhat simpler but closely 
related quantities: the end-to-end correlations of the en- 
ergy density in finite chains, restricting our analysis to 
equal time correlations. 

Since the Hamiltonian involves two kinds of terms, 
J<7q(jI and ha x , to obtain the E-E correlations we 
need to calculate three quantities, < Jo\Oqu\ hi,o\ >, 
< hQO^h L o x L > and < Joi (7 o cr i Jl-\ l^l-i ! >• How- 
ever, the duality transformation simplifies matters, since 
it maps /ioCTq/ilct£ at 8 to J m u^a\ J l-i l^l-i^l at — <^ 
so we only need to compute one of these quantities. Also, 
the mixed correlation function is dual to itself, therefore 
the distribution of < JoiO'o cr i' l i' <7 L > must depend on 
and be the same in the two phases. Note that a 
related single-end quantity, the imaginary time correla- 
tion function < ctq(0)(Tq(t) >, was considered by Igloi, 
Juhasz and Rieger. 14 

The calculation of the quantities of interest requires an 
extension of the methods used so far. The primary rea- 
son for this is that the energy correlations are dominated 
by third order perturbative effects at each stage of the 
decimation, in contrast to the spin correlations, which 
are controlled by second order perturbative effects. 

To be able to carry out higher order RG calculations, 
we develop an approach to the decimation in terms of uni- 
tary transformations; this allows one to follow precisely 
how operators of interest (such as a x ) evolve during the 
RG process. This approach thereby gives a systematic 
way to deal with higher orders perturbative effects even 
in problems previously analyzed using second order per- 
turbation theory. 3 The unitary-transformations method 
is developed in the next section (Sec. II). 



II. UNITARY TRANSFORMATION 
RENORMALIZATION GROUP 

In this section we develop a perturbative scheme based 
on unitary transformations that will allow us to sepa- 
rate the various parts of the Hamiltonian and successively 
simplify the wave functions of the many-body system to 
a hierarchical product of simple spin wave functions. Si- 
multaneously, we must keep track of the original oper- 
ators in order to eventually compute their ground state 
correlations. 

We begin with the first stage of decimation by con- 
structing the eigenfunction of the highest energy part of 
the Hamiltonian and transforming it to take into account 
the low energy parts perturbatively. Specifically, the 
transformation gets rid of the off-diagonal parts that con- 
nect states with large energy differences between them. 
For the Ising chain, this can be done while preserving the 
form of the Hamiltonian. 

Given a Hamiltonian, H, and a many-body ground 
state wave function, \G >, with 

H\G>=E G \G>. (14) 

we can generally make a unitary transformation with a 
hermitian operator S and write: 

e lS He~ lS e lS \G >= e lS \G >, 

(15) 

Heff = e lS He~ lS and e lS \G >= \H > 

with the goal to make \G > close to a product of sim- 
ple wave functions. Such transformations can be used 
to eliminate — or separate — low energy parts in the 
Hamiltonian. 

Let H = Ho + V, where Ho is the high energy part of 
H and V is the remaining low energy parts. The effective 
Hamiltonian is then 

Heff = 7~io 

(16) 

+V + i[S, Ho] + i[S, V] + %[S, [S, Ho]] + 0(S 3 ). 

If we are able to choose S so that 

V + i[S,H o ]=0, (17) 

then H e ff will contain no first order terms. The second 
order corrections to Ho give rise to effective interactions. 
We may now solve for the ground state of H e f / and hence 
the original H: 

H ef f\H >=E G \H > 

(18) 

\G >= e~ %s \H > . 

Iterating this process separates the higher energy parts 
of the Hamiltonian from the lower energy parts. At each 
stage, the effective higher energy parts can be simply 
diagonalized. The remaining non-diagonalized Hamilto- 
nian only has pieces with energy much lower then the gap 
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of the high energy section, which was just diagonalized. 
The ground-state wave function is then constructed per- 
turbatively from the wave- function \H >, which is a hier- 
archical wave-function simply expressible in terms of the 
ground states of the high energy parts of the sequence of 
H e ff's. Each term in the hierarchy will be a spin-cluster 
pointing in the direction of the transverse field. For an 
example see Eq. (35) below. 

Note that this method is related to the flow equation 
approach for interacting quantum problems developed by 
Kchrein and Wegner. 15 ' 16 



A. Unitary RG for Transverse Field Ising Chain 

We now apply the transformations (15) to successively 
reduce the maximum energy scale of the random quan- 
tum Ising Hamiltonian, thereby obtaining a series of low 
energy effective Hamiltonians of the system. We begin by 
choosing the largest energy coefficient in the Hamiltonian 
(1) and denote it Ho (with the corresponding coupling 
the initial energy scale ft/): for example, Ho = —h\of 
(see Fig. 3). Let V designate the part of the Hamiltonian 
that we would like to eliminate. For the above example 
we would like to eliminate the parts connecting site 1 to 
the rest of the chain: 



(19) 



In addition to these two parts, the Hamiltonian also con- 
tains the parts involving the remainder of the chain, 



Hi 



-J-ioV-iOq — ho<Jo —h2&% — J23&2 03 + . • • (20) 



In order to eliminate the first order couplings to spin 
0, S must satisfy equation (17); thus we first choose 



a 2 v 



2/li 



Jl2 
2/li 



<T?<7. 



l u 2 



(21) 



which yields the following terms in the effective Hamil- 
tonian (Eq. 16): 



H 



eff = 

(22) 

Note that site 1 is still coupled to adjacent sites by a 
second order interaction. We would like to restore the 
Hamiltonian to its original form; thus we need to elim- 
inate the new type of interaction. To get rid of it, we 
perform another transformation using 



hoJoi v z 



2h\ ' 



z y 
h°2- 



The effective Hamiltonian now includes 

He// = • • • - J-WV-1&0 — hoCTo 
-J02<J 7 6 (J l (J 2 - h 2<?2 - ^23^1 erf + . . 

-hiof, 



(23) 



(24) 
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H ^e iSl He~ iS i 
eff 



- h -P-l 



-Q 



-w 



J !0^=I„y_ a > J 01 H 1 
h 



o- 



H jj=e iSl He~ iS 2 
eff 



-w 
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FIG. 3: Site decimation. Spin is almost frozen in the x- 
direction due to the strong magnetic field ho. Quantum fluc- 
tuations create a second nearest neighbor effective interaction 
between sites -1 and 1. This interaction is weaker than any 
of J_io, Joi, ho. 



from which we see that in the low-energy subspace of Ho , 
the effective exchange between spins and 2 is given by 



J02 = 



Joi J: 



12 



hi 



(25) 



Since hi is the strongest coupling energy in the chain, 
the resulting effective bond obeys 



J02 "C ho, Joi, J12, 



(26) 



where the sharpness of the inequality is because we as- 
sume strong randomness. We can now partially diago- 
nalize H e ff by writing 



\H>= I -»! |G« > 
\G >= e- iS «e- iS *\H >, 



(27) 



where af \ — >i>= | — >i> and \G^ > involves only the 
spins other than 1 . We are left with a renormalizcd spin- 
chain with the spin at site 1 eliminated, and with an 
effective interaction Jc^Ooci between spin and spin 2. 
[Note that we could also keep the high energy sector that 
involves | <— 1>; the effective Hamiltonian and state of 
the rest of the chain would differ from those of the low 
energy sector because of the presence of in H e ff of 
Eq. 24.] 
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The analog of the above results for the case where an 
exchange interaction, e.g. Ho = — Ji2< J i (J 2^ ^ s eliminated 
is (Fig. 4) 



Ho = - J^oferf 
V = -hiaf - foof 

<? — hl .T y rr z 4- 112 rr z rr v 
t ' a - 2.7,, CT 1 <7 2 + 2.h, a l a 2 



(28) 



2Ji 

S b = -^io 



2 J 2 
ZJ 12 



'1"0 



2Jl2 

op 



'3 U 2- 



This could be obtained by using the duality described in 
the introduction (Sec. IE) and in Ref. 4, or by direct 
computation. The ground state of Ho = — Ji2cfcr| is 
doubly degenerate with spins 1 and 2 either in the state 
I T(i2)>= I Ti> I h> or in the state | | ( i 2 )>= | |i> 
| J,2>- Therefore in the ground state of Ho — of erf spin 
1 and 2 form a ferromagnetic cluster, which we denote 
as (12). We can define cluster operators, and <J* l2 y 
that operate on the spin cluster (12) in the following way: 



Refs. 4,11, with even stronger randomness. The increase 
in the randomness with each step justifies the iterative 
application of the real-space RG as described in this sec- 
tion. In each step we eliminate a high-energy subspace 
of the Hilbert space of the chain, which is gapped by 
2f2 from the remaining subspace. The range of excita- 
tions in the remaining subspace is much smaller than 
Q,. The iterative application of the decimation procedure 
outlined here amounts to separating the Hilbert space 
of the chain into a hierarchy of sequentially decreasing 
energy subspaces. If the coupling distributions expand 



-h a* 



-10-1^0 



-0——0 



H^=e iS2 e iSl He 



- J 0i a a i 



iS i -iS 7 



- J 12°l°2 



'(12) 



^2 



— 0"? <T. 



'(12) 



(29) 



l u 2 



a' 



(12) 



in terms of which 



Keff -H = 
ho(TQ - Jo(12)0O0f 12 ) - h (12) cr ( 1 2) 
- ■7(12)30fl2) a 3 - ^303 



(30) 



with 



h 



(12) 



h\h,2 



(31) 



being the effective transverse field on the new cluster (12) 
that has replaced the pair of spins 1 and 2 that now only 
appear separately in the high energy term in Hq. Again, 
since J 12 is the strongest energy, and strong randomness 
is assumed, the effective transverse field obeys: 



h(\2) < h u h 2 , J12 



(32) 



In both the decimation cases we regain the initial form 
of the Hamiltonian, but with one less spin. As shown in 




- J -J0 a -J a 



FIG. 4: Bond decimation. Sites and 1 are frozen into one 
cluster by the strong Ising interaction, Joi. Quantum fluctu- 



ations produce an effective magnetic field, hoi = 



^01 



which 



flips the composite spin cluster. This field is weaker than any 
of ho, hi, Joi- 



without bounds during the flow, this method is asymp- 
totically exact. 

After applying this set of transformation rules L (the 
original chain length) times, we are left with a single 
spin cluster that carries the moments of some fraction of 
the set of original spins. The ground state of the chain 
is then given by the state in which this cluster points 
in the x direction due to the transverse field. In the 
same way we can also access the various excitations of the 
quantum Ising chain by keeping high energy subspaces in 
the decimation process. 

At the end of the decimation process, the full effective 
Hamiltonian is given by 



Heff = e^e^" 11 . . . e is(2) e is(1) He- is(1) e- isl2) . . . e -^" (33) 



with representing the transformation of the j 'th stage of the original problem by 
of the renormalization. At the final stage the free ground- 
state wave- function, \H >, is related to the ground state .„ m .„ m -aU-D -coo 

|G>=e- lS e- lS ...e- lS e' tS \H > . (34) 
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FIG. 5: Example of hierarchical decimation of a chain with 
four sites. First, site 1 is decimated. At a lower energy scale, 
sites 2 and 3 form a ferromagnetic cluster, which we denote 
(23). Cluster (23) then forms a cluster with site 0. The last 
process is a decimation of the cluster (023). The gound state 
wvae function of this chain is constructed perturbatively from 
the hierarchical wavefunction \H >, which is given in Eq. 
(35). 



Note that the later transformations are the first to op- 
erate on \H >. In the end of the decimation process, 
\H > is a product of cluster wave functions, where each 
cluster points in the direction of the transverse field. For 
example, \H > for a chain of 4 spins as in Fig. 5 is given 

by 



-w 



o==o 



-A 



- J oi°o°i 



- J 12°l°2 




FIG. 6: Evolution of oq in a bond decimation. When the end 
site forms a cluster with its neighbor, the operator ctq gets 
renormalized and gains a factor of hi/Joi- 



\H > = | ^( 23)> I ->1> 

To> I h> I ?3> +| lo> I h> | J.3>) (I Ti> + li>) • in. TRANSVERSE FIELD CORRELATIONS 

(35) 

A. Renormalization of End Spin Operators 



B. Evolution of Effective Operators 

The quantities we are interested in can be written in 
terms of \H > and the set of unitary transformations 
used in the RG process. For example, let us consider 

<o%ot>=<G\o%oi\G>, (36) 

with \G > known in terms of \H >. We can write 

< G\aZai\G >= 

< H\e* s(L> . . . e iS ™e isll) aEaie- isW e- iS ™ . . . e~ isW \H > 

=< H\e^ L) . . . e iS < 2) e ism aZe- is(1) e- is(2) . . . 
e isM . . . e isv e isv axe-isv e -isv . . . e -iS^ \ H > 

= < H\a$al\H > 

(37) 

Thus we see that calculating expectation values of an 
operator A with respect to the ground state \G > is 
equivalent to calculating the expectation value of the ef- 
fective operator A with respect to \H >: 

< G\A\G > = < H\A\H >, 

A = e^ iL) . . . e^ <2) e^ (1) A e^' (1) e"^ <2) . . . e^ s(L) . 

(38) 



In order to obtain the transverse field part of the end- 
to-end energy correlations, < hao^h^a^ >, we must con- 
sider the effects of the two types of renormalization steps 
(decimation of a site or a bond) on the end operators <Tq 
and <r£. In this section we show that the decimation of 
the first bond makes the operator cTq evolve to the ef- 
fective operator (Tg — T^^foiv W e also show that the 

Decimation of the end site yields (Jg = ^t 1 erf , which 
is third order in the off-diagonal terms that couple high 
and low energies in the Hamiltonian. In this derivation 
we neglect all the subleading contributions to the flow of 
the edge operatros; we show in Appendix A that this is 
indeed justified. 

Decimations away from the ends. A renormalization 
step that does not involve the end spin will generally 
leave the end operators unchanged, since the generator of 
the unitary transformation of this RG step, S, commutes 
with (Tq. The exception to this is when the decimation 
involves spins adjacent to the end, for which [S, <Tq] ^ 0; 
however, it can be shown that to all orders in pertur- 
bation theory, there is no contribution to the dominant 
parts of the correlation function from the resulting cor- 
rections to (Tq. 

Decimation of an end bond. If Jqi is decimated, sites 
and 1 will form a cluster (01) (see Fig. 6). The domi- 
nant contribution to the correlation function comes from 
the effective operator cr? 011 . This contribution is obtained 
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from the first order transformation, »S a in Eq. (28): 



0-n = e 



— e iS 



<7n e 



-iS a - 



hi 



(39) 



'(01) 



The first two terms in (39) will not evolve under the 
continuing renormalization. Their expectation values are 



< H\a%\H >=0, 
< H\a$af\H >= 1. 



(40) 



The only piece of (39) relevant to us is the third piece. 
As indicated in Eq. (39), when the operator —<j\<Jq is 
restricted to the low energy subspace in which sites 
and 1 form a cluster, it is equivalent to the cluster oper- 
ator cT(oi) • Therefore when the end spin forms a cluster 
with its neighbor via a bond decimation, the flow of the 
transverse spin is given by 



hi 



-cr] 



J01 



(01)- 



(41) 



End site decimation. If ho is the strongest interaction 
in the chain, site will be decimated. Applying the first 
order transformation S a from Eq. (21) to <Tq yields 



<? — —JQ±rr y rr z 
D ° — 2/ lo °0 Cr l 

~x _ _x I Joi_ n V n z xl _ n x , Jen^z-z 

a a — °b + l l2h a o a n °oJ — °b + h a o a i- 
Applying the second order transformation yields 



(42) 
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[Mo 



im, rr z r 



2h. 



hiJai „V„V , h i J ai 



(43) 



'o u i 



S—< 7 2 ■ 



Once again, the first two terms in Eq. (43) do not con- 
tr 
1. 



tribute to truncated correlation functions since < (Tq > 



< < >- 



0. The flow of (Tq in this case is 



hi J? 



hi 



01 x 

cr-i ■ 



(44) 



eventually combine to form the end-to-end correlation 
function. The evolution of these prefactors is obtained 
from Eqs. (41) and (44) by the method outlined in Ref. 
3. 

At log-energy scale T, we write the effective operator 
of the end spin as 



?5(T)e 



-A(r) 



(45) 



-h a x 



-hrf 




o- 



h j 2 

n i J 01 -A 

? e 
hi 



- J 12°1 Z °2 



FIG. 7: Evolution of ctq i n a s ^ e decimation. Although ctq ob- 
tains an expectation value, its fluctuations are still influenced 
by the state of site 1. This is reflected in the renormaliza- 
tion of ctq as which gives rise factor of Jol a fe ° in the 
correlation functions of a%. 



where <Jq(T) operates on the first spin cluster of the 
renormalized chain at the scale T. Using the logarith- 
mic variables (3 = In ^ , £ = In j , we carL rewrite Eqs. 
(41, 44) and the results described in Figs. (6) and (7) in 
terms of A: 



bond-decimation: A => A + f3\ , 
spin-decimation: A => A + (3i + 2£oi • 



(46) 



The quantities A, (3q, and Iq (length of the cluster con- 
taining the end spin) are correlated at all stages of the 
renormalization. Therefore we must keep track of their 
joint distribution, which we define as 



Evolution of the Correlation Function 



Prob[dlfj,d/3 ,dA] =o;(/3 ,^,A|r)d/3odlSdA. (47) 



As the renormalization proceeds, the effective opera- 
tors <Tq and a x L accrue multiplicative factors that will 



Using the results of ref. 4 and Eqs. (28) and (46), we 
can write the evolution equation for oj((3, I, A|T): 
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du(/3, i, A\r) _ du(f3,l,A\r) 

dr ~ d0 



+ J uj((3 = 0, Z c °, A'|r)P(Co, ll)R{fa.,li)8(\ - A' - 2Co - Pi)6(l -l c -l b - IfiStf - 1 )dl c o dl b dl c 1 dA'd0 1 dC o 
+ J w(J3o, l c , A'\T)R([3 u l<t)P(<; = 0, l b )S(A - A' - (3^6(1 - l c - l b - 1^)6(0 - 0i - 2 )dlZdl b > dl c 1 dA'd0 o d0 1 
- J P(0,r)dl'oj(p,l,A\T). 



(48) 



in terms of the distributions P and R of the log-couplings and lengths of the bonds and spin-clusters at scale L. 

Employing the notation for convolutions introduced in Ref. 4, f(x\) <g> g(x2) = J f(x\)g(x2)S(x — x\ — X2)dx\dx2, 
we can write the above equation in a more compact way: 



du(/3,l,A\r) _ du>(P,l,A\r) 

dr — dp 

+ J lo(j3 = 0, l° c , A'\T) ® P(Co, l b ) 4 R{/3, lf)5(A - A' - 2<„ - P)dA'd( 

+ J w((3o, l c , A'\T) ® R(0 1 ,l c 1 ) ® P(Co, l b )S(A -A'- fa)dH 
-J P(0,r)dl'u>(f3,l,A\T). 

I 



(49) 



By Laplace transforming oj(0, I, A) with respect to 
both I and A, 



uj(0,y,X) = JdlJ 
we obtain 

du(f3,y,\\r) 

dr 



dAe- ly e- XA Lu{f3,l,A) 



(50) 



dp 



+cj(0 o = 0, y, X\T) J e-^ x P(C, y)dQR{0, y) 

+ / y, X) I R(0 1 ,y)P(Co = 0, y) 
-P(O,y = O)u;(0,y,X). 



(51) 



In Ref. 4 the scaling limits of the functions 
R(0, y), P(/3, y) are derived: 

p(C,y|r) = T( y ,r) e -^ r ), 

R{P,y\T)=T(y,T)e-Wy?\ 

T(y,r) 



Afa) r sr 



T(y,r) 



ahA(y)r l 



-5r 



sinhA^r'' ' 

T(i/,r) = 5 + A(y)coth(A(y)r), 
u(y,T) = -8 + A(y) coth (A(y)L) ; 



(52) 



My) = Vv + ^- 

Using these results and the corresponding notations, we 
get 



dr ~ ' ap 1 '"V"' 2c+m 
+T(y, r) Ju((3 ,y, X)e- T ^-^T(y, T) 
■5(/3 - O ~ PiWodPi 
-T(0, !>(/?, y, A). 



(53) 



To solve this we write co(0, y, A) in the following form: 

co(p, yi X) = W(y,X)e-^ T , (54) 
under which (53) becomes 

^ = -(r + X)W + ^^lw-mT)W. (55) 

(XL — 1~ U 

Using the definitions of the functions T, T, u, and r , we 
can integrate Eq. (55) and find 



sinh(A(y)rj) / 2\+u(yX I ) \ r(0J) 
' sinh(A(y)r) \ 2X+u(y,r) J r(0,r 7 ) " 



(56) 



The value of Wq can be found from the normalization 
condition on the distribution co(0,l,A). To be a properly 
normalized, u>(0,l,A) has to obey 



d0w(0,y = 0, A = 0) = 1. 



(57) 



This implies W = r(0,r 7 ). 

Before proceeding, we should check what initial con- 
ditions this distribution satisfies. Setting r => Tj and 
y ^> we see: 

u(0, 0, A|rj) - t(0, r J )e- r ^' r ^ /3 e- A/3 (58) 

The inverse Laplace Transform in A yields 

u(J3, 0, Air,) = r(0, r J )e-^°' r ^^(A - 0). (59) 

The delta function in Eq. (59) shows that the initial value 
of the correlation variable is the same as the transverse 
field, as it should be for this part of the energy - energy 
correlations. 
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C. Last Decimation Step 

After carrying out the decimation process L — 2 times 
and following the flow of the edge energy operators, we 
end up with two clusters. The remaining clusters corre- 
spond to the left and right edges of the chain, and each 
has a distribution u(/3, l,A\T) associated with it. In the 
next decimation step one of these clusters is decimated, 
and a single cluster forms. At this stage we can compute 
the distribution of the energy correlations from the flow 
of edge operators and the remaining couplings. Note that 
the r at which a single cluster forms is logarithm of the 
gap between the first and second excited states (as was 
also noted in Rcf. 14). 

The computation of the truncated correlation function 
is as follows: 

C X L X =< G\h a x h L a x L \G > 
- < G\h a x \G >< G\h L a x L \G > 
= < H\e- Kt afafe- K ^\H > 
- < H\e~ Ae af\H >< H\a x f e-^\H >, 

where Kg and A r are the correlation factors picked up in 
the RG process, Eq. (46), for the left end and right end 
transverse spins respectively, and we have labeled the last 
remaining spin clusters i and r. The correlation function 
can also be written as a sum over excited states: 

Cl x =< h a x h L al > - < h a x >< h L a x L > 
= < G\(h a x - < h a x >){h L al- < h L a x L >)\G > 
= £ <G\hoo$\il>>< ^\h L cjl\G> 

= e -A e -A r £ < H\a x \ip >< ip\a x \H >, 

(61) 



where the sum over ip runs over all states except the 
ground state. 

At the penultimate decimation step two processes are 
possible: 

Bond Decimation. The two remaining clusters be- 
come a single cluster, and the only remaining coupling 
is the transverse field, which makes the combined cluster 
point in the x direction: 

W (L) >= -^(l T*> I Tf> +1 li> I ir>), (62) 



where I and f represent the end clusters at this final stage 
of the RG. This yields (60) 

C xx = e -A e -A r = e -A. ( 63 ) 

hence 

A = A r + A^, (64) 

where A is the desired log of the energy-correlations. 

Site Decimation. In the case of a site decimation it 
is unimportant which of the last surviving clusters gets 
decimated. Let us assume that it is the left cluster, I. 
This involves the unitary transformation (dropping the 
tildes) S a = ^^vj^h an( ^ ma kes the ground state \H > 
be 

n {L) >= | ->/> | ^ r >= u> +1 U>)(\ Tr> +| lr>). 

(65) 

By using the sum form in Eq. (61) and applying S a we 
get 



C xx e Al+Ar w < H\e lS «a x e- lS «\ip >< ^\e lSa a x e' lSa \H > 

(66) 

i>^G 1 

This yields 

A = A r + A/ + 2C/ f . (67) 
The analog of Eq. (25) in Ref. 3 for the measures defined here is 
rfProb(A,, l t , Pi, lid, 11, 01, ll C2, • • • ,lr, fa, K\L, T) = 

aru(Pe, At, k)P{<;ul b 1 )R((3i,l c 1 )P{(2, l b 2 ) ■ ■ ■ w(/3 r , A r , l r )8(h + l\ + lt + l\ + . . . + l r - L)d{(3 i }d{C i }dK r dK e .d{l i }, 

(68) 
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with: 

oo 

= J( p (0,l) + R(0,l))dl. (69) 
o 

Let us now define the function J(A, T|L) as 

dProb[chain of length L becomes a single cluster at T with log(Cx) = A] = J(A, T\L)dTdA. (70) 
The probability distribution of log(Cx) is then given by 

oo 

f(A,L) = J J{K,V\L)dT. (71) 
o 

The function J(A, T\L) has two contributions, the first contribution cames from the case of the penultimate dec- 
imation being a bond decimation (Eq. 64). The second contribution comes from the case of a site-decimation (Eq. 
67) . The combination of the two contribution yields: 

j(A,r|L) = 

a r J P(0, l b er )u)((3 e , At, k)tj{(3 r ,A r , l r )S(k + l\ r + l r - L)S(A -A t - A r )dAedA r dledl b er .dl r d/3 r d(3 e (72) 
+2a r / P(( r , l b £r )u(0, At, lt)u(0 r , A r , l r )5(k + l\ r + l r - L)5(A -A e -A r - 2C r )dAtdA r dkdl b lr dl r df3 r dC,r 

where l b r is the length of the effective bond connecting the last two clusters. 
The Laplace transform of J is considerably simpler to write: 



J(A,r,y) = 



a r 



P(0, y) ( J u{(3, X, y)d0) 2 + 2 J A P(C, y)d(uj(0, A e , h) J uifi, A, y)d(3 
fl r^(0, y) (tShIa) ((r(»|r)+A) + 2 («( a |r)+2A)) ■ 



(73) 



Eq. (73) is one of the main results of this paper; from 
it we will derive the typical and average correlation func- 
tions, as well as information on the distribution of Cl- 



IV. EXCHANGE ENERGY AND CROSS 
CORRELATIONS 

A. Boundaries and Duality 

The calculation of the hal^ — h<j rj [ correlations in the 
previous section simplified greatly because the operators 
(Tq l considered were end operators. In this section we 
calculate expressions for the end-to-end correlations of 
the exchange energy density Jo z o z . We consider the 
special case for which one or both end-transverse-ficlds, 
ho and hh, are zero, which simplifies the calculation con- 
siderably. We will argue that the universal features of 
the correlations will be the same as in the general case 
with non-zero end-fields. 

The simplifications with vanishing end transverse- 
fields arise because this makes the exchange energy be 
an edge operator in the sense that it is the first and last 



energy operator in the Hamiltonian: 

H = — JoiUquI — h\o\ — .. .-h L -\o x L _ x - Jl-il<7l-i<7l- 

(74) 

And edge bond is the dual of an edge site, and there- 
fore the edge-to-edge exchange- energy correlations of the 
ground state of the Hamiltonian in Eq. (74) are duals 
to the transverse-spin correlations calculated in Sec. III. 
This is explained below. 

As was explained in the introduction (IE), the duality 
transforms a bond to a spin and vice versa. In the pre- 
vious section we considered a chain that terminates with 
a site that has a finite transverse field on it, ho > 0. The 
dual of this edge is a chain edge that terminates with 
a nonzero bond, Jq'v = h . The site 1' is the dual of 
the bond J i and therefore experiences a field hy = Joi- 
Since there is no bond J_i , i.e., J_i o = 0, the field on 
site 0' is zero as well (see Fig. 2). 

In what follows we calculate the end-to-end correla- 
tions of the exchange energy and the cross-correlations 
between the transverse spin and the exchange energy. In 
both cases we will assume that the chain terminates with 
the energy operators whose correlations we calculate (as 
in Eq. (74) for the exchange-energy correlations). In the 
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case of exchange energy correlations both edges of the 
chain we consider will terminate with a vanishing trans- 
verse field, ho = h.L = 0. Similarly, when we calculate the 
edge correlations between the exchange energy of sites 
and 1, and the transverse spin on site L, the edge trans- 
verse field ho is set to 0. These rules allow us to use 
the dual of the function uj((3,1,A\T) which was derived 
in Sec. IIIB. We define the function <f>((,l,A\T) as the 
dual of 1, A|T). This function will keep track of the 
correlations and evolution of the operator Joi< J o cr i m the 
same way that uj((3, I, A\T) was used to keep track of the 
correlations and evolution of the operator hoa^ . 

The calculations carried out in this section assume that 
one or both edge transverse fields are zero, but this does 
not limit the generality of our results for universal quan- 
tities. We expect that when the edge transverse-fields 
are non-zero, the correlations of the last bond in a chain 
will only be modified by a non-universal multiplicative 
factor from the correlations in the special case with no 
transverse-field on the end spin. 



B. Evolution of Edge Exchange-Energy Operator 

As discussed above, the edge exchange- energy Joid^af 
is dual to the and edge transverse field operator hoa^. 
Therefore we can obtain the distribution function for the 
evolution of the edge exchange energy operator from the 
results of Sec. III. 

By making use of the duality (Sec. IE), we can trans- 
form all the results obtained in Sec. Ill to the dual chain. 
As stated above, we define the analog of w(/3, I, A\T) to be 
<f)((, I, A|T): (/>(£, I, A|T) keeps track of the bond strength 
of the end bond, its length (including the length of the 
h = end site), and the log contribution to the correla- 
tion, A. <^(c^, /, A|r) is obtained from the dual of (56): 

<KC,y,A) = $(y,A)-e-< A -^> r ) 



A) = e<A-*>Cr-r,) 3gg$«(0, r). 



(75) 



C. Exchange Energy Correlations 

The results for the exchange energy correlations are 
given, by duality, by Eq. (73), with 5 — > —S. This yields 



J(A,T,y) 



a r^(°, V) (uSrH-A) ((«(y|r)+A) + 2 (r( 1 ,|r 1 )+2A)) 



(76) 



D. Cross-Correlations — Last Decimation Step 
and Final Expression 

In order to obtain the cross correlations we need to 
combine the results for the edge transverse-spin flow and 
exchange-energy flow. In analogy to Sec. Ill, putting 
together the two flows happens in the penultinate step 
of the RG flow. The accumulated multiplicative factors, 
along with the couplings of the renormalized chain just 
before it is completely decimated, will determine the to- 
tal correlations between the transverse-spin and exchange 
energy. 

In contrast to Sec. Ill, the last needed step of the RG 
to obtain the cross correlations involves the transverse 
spin of one of the two clusters, and the bond between 
them: 



C x-B 



L -<G\hoolJ L - 1L oi_ 1 oi\G> 
< G\h a%\G X G\J L ^ L al_ x al\G > 



-< H\e- Ae a^a z ,aie- A -\H > 



(77) 



e r 

< H\e- Ae af\H >< H\a\a z f e 



\H>, 



where B stands for bond. The two possibilities for the 
last step of the decimation process are the decimation 
of the bond (Jji f ), or of the t cluster (hi). These two 
processes are dual to each other; hence we only need to 
consider one of them. Let us consider the site decimation. 
As before, the ground state will be 

\H > = | | ^r>= \{\ Tl> +1 ll>)(| Tf> +1 lf>). 

(78) 

From the transformation S a = ^^' G f IJ \ that induces 
this decimation, the correlations are found to be: 

C x-B e A e +A r =< H\e iS -afe- iS -e iS -afa z f e- iS -\H > 

- < H\e iS «o^e- iS °<\H >< H\e lS *<jf(jfe- lS «\H > 

=< H M - t°l°?) W + t°i )lH > 



■1m 

hp 



-e 



(79) 



Thus the cluster decimation process yields 

A = At + A r + Q f . (80) 

By duality, the bond decimation process yields 

A = A e + A r + (3- t (81) 

Following the reasoning that led to Eqs. (72, 73) we 
get: 



Since the results for the exchange-energy correlations are 
identical to that of the transverse correlations, we will 
only analyze the later. 



J(A,r|y) = a r w(A,rXA,r) 



1 



+ 



r(y|r) + A u(y\T) + \J 
(82) 
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This is the second main result in this paper, and it is 
analogous to Eq. (73). Here we must bear in mind that 
the correlations obtained here are negative (see Eq. 79). 
This is to be expected, since the two operators, a x and 
<j z a z 7 try to impose competing orders; One tends to dis- 
order the system and the other to order it. 

V. RESULTS 

The above results (73, 82) in principal allow the calcu- 
lation of the distribution function for the log-correlations, 
A = — log(Ci), of long finite chains. In the following sec- 



tions we calculate the average correlations, Cx, and the 
distribution, f(A\L), for all 5. 

A. Average ha x - ha x and Jo z a z - Jo z o z 
Correlations 

Derivation of the Average. In this section we derive 
the average x-x correlations. The BB correlations of the 
exchange energy are obtained from the x-x correlations 
upon the transformation S — ► —S. 

In order to obtain Cff, we begin with Eq. (73) in the 
following form: 



J(X T\v) = T( "' F) sinh2 < Ar j) A ^A-a+AcothtArj)) 2 (2A+3<5)r+2(A+5)r r 
v ' I"' u(Q,T) sinh(Ar) (26— X) 2 



(2(5— A) 2 



1 . 2(25- A) sinh(Ar) 1 

P+A) sinh(Ar)+A cosh(Ar)] 2 + [(2A-<5) sinh(Ar)+A cosh(Ar)] 3 [(2A+<5) sinh(Ar)+A cosh(Ar)] 2 

First, we perform an inverse Laplace transform in y to recover the length dependence: 

J(\T\L) = £ ( _i) ^(o,r) ((2A-J)r f +i)' e - ( 2A+M)r+2(A + i)r J e -(V+(^)> 



4L(ri7r) 4 
K 

2(»i) ! l 



u(0,r) (2(5— A) 2 

n— 1 

P (A + 5)r;i i 8L (ktt) (25-A) _ 
e -I- r 9 (2A-5) 3 6 



4L(«7r) 4 



r6(2A-5)(l+ rJ(2A - 5) 



IL__ 

2(n7r) 2 £ 



_ e (2A-a)r 3 



(83) 



(84) 



This is obtained by approximating the roots of 

a sinh(Ar) + A cosh( Ar) (85) 

with a w 1, by y„ = — 5 2 — (-^j?-) 2 (l — j^) and expanding 
(85) around these roots: 

asinh(Ar) + Acosh(Ar) w V(-l)" +1 »(y - Wn )^. 

n=l 

(86) 

The roots in Eq. (86) are given as an expansion in powers 
of 1/r. Since we are interested in T > 1 we are content 
with only the first two terms; in fact, as can be seen by 
the following Eq. (87), only the first nonvanishing power 
of 1/r contributes to the average correlations. In addi- 
tion, Eq. (86) is only valid for y n <C 1, i.e., for nn < T. 
But since we are interested in the large length behavior of 



the correlations, we can restrict our calculation to small 
values of n, as they give the slowest decaying term in the 
correlations. 

The desired result is obtained by performing the T in- 
tegral. This integral is dominated by the large exponent 
in 

e -(i ¥ L) 2 L _ ( 2A+,5)r(+o(iog(r))+o(Lr- 3 )) = g g(r) ^ 



which has a saddle point at Ts = 




The 



exponential dependence then becomes 

e9 (r )Rje ^)- 3 (^) 1/3 ( r -^ 
and the saddle point integration yields 



J(A, L) M £ (-l)" +1 2 2 / 3 yf (^)l/3_l_ e -5 2 L-3L 1 / 3 (^) 2 / 3 (A+5/2) 2 / 3 



71=1 

vr 2T,(\+&) ( 2Arj+A 2 i (89) 

Ae V. x - 25 J (2A+5) 2 / 3 

X (c 2 -^ {2X+5)3 I c l-^x ^-A)(2A+j) 2 (2X+5? 
A I e (A+<5)(3A+2<5) ^ L (2A-5) 2 e (4A)(2A-<5) 
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This result is valid for the critical regime (4? ^> L) when 
jj/3 <A<I, and away from the critical regime <C 
L) away from A = 0. To get the equivalent expression 
for A => 0, we need to be more careful with the third 
term of equation (89) and get the next order corrections. 
Note that there is no singularity in this expression at 
A = 25; this will have implications for the off-critical 
large L behavior. 

Result. To get the final result for C£ x all that remains 
is to set A => 1. Neglecting terms suppressed by factors 
of 8 <C 1 we obtain 

Cf K ^x_l_ e - 4 ^-3 i V. W V3 (1+4/2) VS i (9Q) 

Also, for the exchange-energy correlations we obtain 

A** _L e -« a i-3LV. W Vs (1 _ 4/2) V. ) (91) 



Result. The average x-B correlation is obtained from 
the above by setting A = 1: 

Cf* » -Al-B^e-* 2 ^ 1 ^ 2 ' 3 , (94) 

which is almost the same as Eq. (90) but the 8 depen- 
dence of the above result is strictly symmetric with re- 
spect to 8, as expected for an object that is self-dual. 

C. Typical Correlations 

One of the striking features of random quantum sys- 
tems is that typical correlations are usually very different 
from average correlations. Average correlation functions 
can be, as here, dominated by samples (or spatial re- 
gions) with anomalously strong correlations. The typical 
correlations are much smaller and decay faster with dis- 
tance. Indeed, for the random Ising chain, the typical 
correlations hold for almost all long-but-finite samples. 



where A xx and Aq z are non-universal coefficients. 

When the chain is not critical, we notice that the expo- 
nential decay is controlled by the same correlation length, 
£ = 4r, as the order parameter correlation function. 



B. Average Cross Correlations 



Derivation of the Average. In complete analogy with 
the derivation of the previous section, we proceed from 
Eq. (82) in the following form: 



(92) 



(93) 



The average end-to-end correlations are dominated by 
extremely rare samples. 

In the off-critical regime the typical correlations are 
well characterized by the average log-correlations 



Ctypical e — g 



However, at the critical point — log(C) = —A is of the 
same order as the logarithm of the typical correlations, 
but the typical correlations will also have a very wide 
spread. More precisely, there is a random proportion- 
ality Constant relating the log of the correlations to its 
average; this constant is random and widely varying. We 
will first investigate the typical correlations at the critical 
point, and then consider the off critical regime. 

Typical Correlations at the critical point. The average 
log-correlations A are easily found at the critical point. 
Going back to Eq. (73) and setting 8 = 0, we see that 
the Laplace transform in L of the x-x log-correlations is 
given by 



Li" l " 



J(X,r\y) « (fi^±i^n e - 2 .(r-r,)^ 
x ( _L i l _± I 

\(A+25) ((A-5)sinh(Ar)+Acosh(Ar)) ^ (A-25) ((A+5) sinh(Ar)+A cosh(Ar)) 

I (A-g) 1 , (A+5) l 

' 5(25— A) ((2A-5) sinh(Ar)+Acosh(Ar)) ^ 5(25+A) ((2A+5) sinh(Ar)+A cosh(Ar)) / 

Performing an inverse Laplace transform in y and then performing a saddle point integration in T, we find: 
J(A|L)« £ e ^ J ((2Ar J+ ig-(5r f )-) 2 ^ ( ^ )1/ 3_i 7 _ e -5^-3^/3 ( „ 7r) V3 A2 /3 



71=1 



; + e 2 ~ A+? _ (A-5)e + 2A^7 (A+5)e 

A-25 + A+25 5(A-25) + 5(A+25) 

This result, as well as Eq. (89), is valid for jjjj <C A <C L for all small or zero 8. 
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7 2tanh(r v ^)(3(r-3r J ) v ^+10tanh(r v ^)) , p 



1'j 



(95) 



2 

,,3/2 



^ S cosh 4 (a:) 



dx 



7tt 

iy 3/2 



By performing the inverse Laplace transform we get 

L 4 W (96) 

w 3.1VI. 

The result in Eq. (96) should be compared with the 
critical behavior of the average correlation: 

log (~) « -37r 2 /3LV3. (97) 

The typical correlations decay as ~ e~ fcv ^ with fc of order 
unity but random with a computable universal distribu- 
tion. But the average correlation function is <~ er c L 1 , 
which decays much more slowly. As claimed above, this 
means that realizations of the quenched randomness that 
have an exponentially low probability dominate the av- 
erage. 

By differentiating (95) once more with respect to A, wc 

get (A**) . Using that we get for the standard deviation 
of K xx : 

\J {A xx f - w 5.6\/I. (98) 



The distribution of A/y/L is thus indeed non-trivial for 
long critical chains. 

For the cross correlation function, Ja z a z — ha x , the re- 
sult we get for the average log-correlation is (by a similar 
calculation) 



A r s «37?vi+o(i) 

« 3.0VX, (99) 
\] (A x ~ B f - ~A x ~b 2 S3 5.4VX. 

Note the similar — but not identical — behavior of the 
two results (99) and (96). 

Off Critical x-x Correlations. To investigate the x-x 
energy correlations in the off-critical regime, we pursue 
a different course of action. Instead of setting A to 1, wc 
invert the Laplace transform with respect to A in expres- 
sion (73) and obtain J(A, T\y) in terms of A. 

Equation (73) can be written in the following form: 



J(X,T\y) = ggff ffiffy A(2A-5 + Acoth(A r/ ))2 e -2A(r-r,) e+2 ^- 3 ,r 

1 | 2 

(35+ A coth(Ar)) 2 (5+A+A coth(Ar)) 2 ~T (3.5+A coth(Ar)) 3 (.5+A+A coth(Ar)) 

4 4 

(3.5+A coth(Ar))(2A-<5+A coth(Ar)) 3 (3.5+A coth(Ar)) a (2A-<5+A coth(Ar)) 



(100) 



This form of J(X, T\y) lends itself to inverting the Laplace transform and recovering the A dependence. This gives 
(neglecting Ti, as before) 



r(0,r) 1 
u(0,r) sinh(Ar) 3 



A 3 e+ 2 ^- 35r e(A-2(L-F / )) 



(A-2(r-rj)) , 2 \ _-(A-2(r-r 7 ))(5+Acoth(Ar)) 

^ (3.5+A coth(Ar)) 2 T (3.5+A coth(Ar)) a J c 

i(A-2(r-rj))(-5+Acoth(Ar)) 



(101) 



, I 2(A-2(r-r 7 )) 2 4 

" r ^ (3.5+A coth(Ar)) (35+Acoth(Ar)) a 



e 2 
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with the Heaviside step function. Off critical, for long enough chains, specifically with, L 3> £ ~ ji and the 
concomitant log-energy scale T8^> I, one can expand 



5 + 



26 



(102) 



and hence obtain 



J(A,r|j,)«|<J| 3 e+ OT '-(WI> r e (A _ 2r ) 



(A-2(r-r I )) 

(45? 



/ 2(A-2(r-r I )) 2 
+ I m 



-(A-2(r-r 7 ))(«+|«|) - 5 ^ T (A-2(r-r J )) 



(103) 



e 2 



|(A-2(r-r I ))(-5+|5|), 



(A-2(r-r 7 )) 



From the simple form of the y dependence in this limit, one can invert the Laplace transform by inspection to obtain 
the L dependence. 

In the paramagnetic phase (S > 0), J(A, is sharply peaked for long chains at 2(r — Tj) + 45L. In the 
ferromagnetic phase (6 < 0), it is instead sharply peaked at at 2(r — Tj) + 2|d|X. Integrating over T gives exponential 
decay in L. Thus the distributions of the end-to-end transverse field log-correlations of long off-critical samples will 
have the form : 



oo 

f xx (A,L) = J J(A,T\L)dT 



e -25(A-45L) e ( A _ 4(5i ) 6>Q 

e -l«l(A-2|«|i,) e ( A _2| ( y|£) A<0 



(104) 



[Recall that rj is the Heaviside step function.] This behavior, with exponential decay of almost all samples with 
a characteristic length that is much shorter than the correlation length is similar to that of the order parameter 
correlations in the paramagnetic phase as discussed in Sec. I C of the Introduction. 



Off Critical x-B Correlations. The same analysis can 
be applied to the x-B correlation function. The math- 
ematical expressions are simpler, but the result is more 
interesting. Since this correlation function is symmetric 
with respect to 5, we can choose to carry out the anal- 
ysis in the paramagnetic phase, 5 > 0. Using the same 
simplifying limit as before (Tj 0, TS 3> 1, y <C <5 2 ), 
and keeping only the dominating terms in the disordered 
phase, Eq. (82) becomes 



J(X,T\ y )^e- 2X ^- r ^- 2Sr 5 2 



Vv 
\ y 

A ~ 45 



Viz 

A ~ 25 



(105) 



The inverse Laplace transform in A and y of this leads to 



(0[A-2(r- 



Fj) - 26L] - 



-25rr2 



e[A-2(r-Fj) 



A5L]). 
(106) 

By integrating over T we get for the distribution of the 
end-to-end log cross-correlations, A = A°^ B , of long off- 
critical chains, 



oo 

f(A,L) = J J(A,T\L)dT • 



(0 A<2\S\L 

1 _ e -5(A-2|6|L) 2\S\L < A < i\S\L . 

{ (l_ e -25 2 £) e -|5|(A-4|5|L) A>4|«5|L 

I 



(107) 



Most long samples will have 26 L < A < A8L\ remem- 
bering that this expression is valid only for Ld 2 ^> 1, 
we see that the distribution will be roughly constant in 
this range. For larger A, the distribution decays expo- 
nentially. 



In the ferromagnetic phase, the same result for the 
cross correlations will obtain with S replaced by \5\. The 
behavior in this phase contrasts with that of the x-x cor- 
relations whose distribution of A is peaked near 2\5\L and 
thus are typically stronger than the cross correlations. 
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D. Energy - Energy Correlations and the Energy 
Gap 

Looking at the earlier results for the average of the 
energy gap, AE, of finite chains, 3 we observe a strong 
resemblance to the results obtained here for the average 
E - E correlations. In particular Eq. (60) in Ref. [3], 

AE ~ L 1/6 e~H~J , (108) 

gives the average gap at the critical point S = 0. We see 
that 

2 log AE ss log Cp (109) 

Some relation between the gap and the energy correla- 
tions is to be expected, but the behavior of the two quan- 
tities is surprising in the degree of the similarity. We will 
see that this relationship between the gap and the energy 
correlations arises in the structure of the RG flow. 

In Eq. (103), the Heaviside function 6 (A - 2(r - Tj)) 
implies that — log(C^) > 2(r — Fj). Since the appro- 
priate T at which the decimation establishing this corre- 
lation occurs is T = — log AE, this implies that 



L 




(110) 



To understand this inequality, consider the correlation 
coefficient e~ A associated with an end site, and compare 
this to the transverse field h on the end spin cluster. The 
strongest correlations will occur if the chain undergoes a 
series of bond decimations. Looking at Fig. (6) one can 
see that in this case the evolution of e~ A and h are exactly 
the same; in each end-bond decimation they acquire a 
factor of the at that scale. The strongest correlations 
dominate the average correlations. Therefore chains in 
which the energy-correlations and the gap are strongly 
correlated also dominate the average correlations. 

VI. CONCLUSIONS 

In this paper we have investigated the various con- 
tributions to the end-to-end energy correlations of ran- 
dom quantum Ising chains in the universal regime of long 
chains near the quantum phase transition. In principal, 
the main result obtained here, the Laplace transform of 
the distribution of the logarithm of the correlation func- 
tions (73, 82), can be used to obtain the complete distri- 
butions in the scaling limit. We have explicitly computed 
the average and typical correlations, as well as some other 
aspects of the distributions, in various limits. The aver- 
age correlations are dominated by exponentially (in the 
chain length) rare samples. Nevertheless, they still de- 
cay as e~ CL 1 at the critical point. This is in contrast to 
the power law decay of the average order parameter cor- 
relations. The various components of the energy corre- 
lations are qualitatively similar, although their distribu- 
tions differ. The cross correlations between the ordering 



operator, Jo~ z o~ z , and disordering transverse-field opera- 
tor, ho~ x , are negative because of their competing effects; 
they are also self-dual. 

The average correlations in the off-critical regime decay 
with the same correlation length: £ ~ ^ , as the order pa- 
rameter correlations. The typical correlations, however, 
decay much faster, and have the same correlation-length 
exponent, v = 1, as the pure system. In Sec. (VC) 
it was shown that the distribution of Cgcr^ is strongly 
peaked near exp(— A5L) for 6 > and near exp(— 2|5|L) 
for S < 0, indicating a surprising asymmetry between 
the two phases. This asymmetry is a result of the differ- 
ence between the ordering and disordering components of 
the energy density and their behavior in the correspond- 
ing phases. At the critical point, the typical correlations 
decay as ~ e~ KL 1 with K a random variable; this is 
similar to the typical order parameter correlations. 

The behavior of the end-to-end energy correlations 
turn out to be similar to that of the energy gap. We ex- 
plain this in terms of the rare realizations of the quenched 
randomness that dominate the averages: these are such 
that the gap and the correlation function involve essen- 
tially the same product of ratios of J's to ft's. 

Unlike previous RG calculations of properties of the 
random transverse-field Ising model, the energy corre- 
lations required the development of a formalism that 
goes beyond second order perturbation theory. Perform- 
ing the RG transformation by unitary transformations 
proved to be a useful tool that allows one to follow read- 
ily the evolution of effective operators. Here we have 
focussed on end-to-end correlations because these are far 
simpler to handle analytically: correlations in the bulk 
of the chain involve effective operators on both sides of 
the objects of interest and are much harder to deal with. 
Nevertheless, they could be computed by keeping track 
of the needed distributions numerically. 

The unitary transformation RG method can also be ap- 
plied to the calculation of correlations in imaginary time. 
In this case, primarily average quantities have been cal- 
culated thus far (for instance, see Refs. [14,17,18]) but 
progress on distributions should be possible utilizing the 
procedure described here (Sec. II). The present method 
may also be applied to other models both in one dimen- 
sion, and, numerically, in higher dimensions. 
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APPENDIX A: EFFECTS OF END OPERATORS 

When calculating the energy correlation functions we 
kept terms that were third order or higher in the per- 
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turbation expansion. A conccern to the validity of our 
treatment is that terms that were produced as fourth or 
fifth order terms at an early stage of the RG flow, T\, 
become more relevant than terms that we kept that are 
produced at a later stage of the RG flow, T 2 > IV In 
this Appendix we verify that end operator effects that we 
excluded in the text can not give rise to leading order con- 
tributions to the various computed energy correlations. 

In addition, The flow of the operator ctq when the end 
spin is decimated, will include a term a§af . Naively, as 
this is an order — rather than a disorder — energy op- 
erator, it may make the disordered and ordered phases 
look indistinguishable as far as the transverse spin cor- 
relations arc concerned. Wc will show here that this is 



not the case, and although these extraneous operators do 
appear, their contribution is at best subdominant. 



1. End Site Decimation 

In order to prove the above claims, we need to investi- 
gate the additional operators that arise when decimating 
ends. Revisiting the process of decimating an end-site, 
we consider the flow of the operator hoa^. The follow- 
ing table describes the series of transformation, and the 
effective operators that contribute to h <jQ-. 



S a : 
S b : 
S c : 
S d : 



2h "0 °1 

27^— a a l 2 

JoiJl2hih,2 z x —V 



Cr (T l Cr 2 



-Joicrgcrf 
Jaih 1 y „y , Jmh 



Jo 1 Jl 2^-1 



2— Ci 



(Al) 



'0"1 u 2 



The fourth line and the second line contain the terms 
that are most likely to give a leading contribution to the 
correlations; these operators have a non- vanishing expec- 
tation value in the ground state of the decimated part. 
Keeping them all we have 



T 2 h T 2 

J 01 n l x | % t .2 zz 

1 + TA 12 1 1 2 
n 



I 2 
J oi 



hi 



h, 



4 J\2h\h2u\a\. 
o 



(A2) 

Another process that we need to consider is a bond 
decimation close to the end. If we decimate the first 
bond, — Joicr^l, then /iqCq becomes 



h a <7% 



hght 



(A3) 



-J\2hi(jlala\ = 

/l (oi) cr foi) + w h {oi) h 2J(ai)20- v {m) o-l, 

where the parentheses signify effective spin clusters. Go- 
ing further away from the end, a decimation of the second 
site in the chain, —h\af, will give 



hoo-Q + j?Jo2hoh 2 o-o<T%, 



(A4) 



with site 1 eliminated. Note that site is still the first 
site but site 2 is now the next. 

From the above processes a pattern emerges. The fifth- 
order perturbation calculation above produces two dan- 
gerous operators: 

a X, ola\. (A5) 



We need to verify that these operators do not pro- 
duce leading-order correlations. The first step is 
to observe that instead of seeing the bare opera- 
tors, Ctq, o-Qcrf, o~qO~\, appearing with varying pref- 
actors, we see them appearing in the combinations 
ho&Q, ■j^2h 2 ) Joi<TQ<jf, ■^ 2 -hoJoihio~QO~\, where fl is the en- 
ergy scale at which the operator appeared. There may be 
additional prefactors, which we will consider shortly, but 
first let us establish that these forms have an invariant 
structure. 



2. Invariant Operators 

In this section we will show that the form of the oper- 
ators defined as follows: 



hoo-Q, 
h Joihi<jQ<j\ 



(A6) 



is preseved during the RG flow. The coefficients ho, 
/iqJoi, and hoJoihi will only be replaced by their renor- 
maized counterparts every time that a decimation affects 
them. This is demonstrated in teh folloing examples. 

The simplest example is the operator /loCo- Under the 
decimation of — JoictqCtJ , 
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hghl y y _ ,eff X 

n<oo- => — <r a 1 - h, 01) a, 01 y 



Jo 



(A7) 



Another example is the operator ^% 3q\Oqo\ . The decimation of the second bond, — J^ofcrf, modifies this operator 
as follows: 



In the case of the operator fp-ft-o^oi^-iCgCrf, the corresponding transformation S a = jj^^Wt yields 

haJoih^lal => h J 0{12) ^j^-af,af a\ = h Jo(i2)^(i2)^o cr ('i2)- 



(A8) 



(A9) 



Decimating the second site in the chain, —h\u\ involves the transformation S a = xf-of erf , yielding the following 
flow for the two operators J<no%o\ , jp-fto Jmh\a v ^o\: 



l2 t zz . h 2 ^01^12 z x z 



i 2 ref f z z 
T ^O' 7 !^ = /l J 02 ^O^' 



hi 

The second transformation in this same decimation process is Sb = Jl ^i 2 afa 2l which gives rise to: 

h a J 01 hia^af =4> h J °] Jl2 hi^a%afa% = h^J^ h 2 ala\ . 
hi tii 



(A10) 



(All) 



r 



In all cases the flow due to decimations leaves the three 
forms of the end operators invariant. We excluded here 
the cases of a decimation of the first site or bonds; these 
are considered below. 



3. Displacement Prefactors 

As mentioned before and seen from the results of (Al), 
there are still multiplicative prefactors coming before the 
invariant operator forms. In Eq. (A2), for instance, all 

three operators from Eq. (A6) have the prefactor (7^) • 

This suppression can be associated with the displacement 
of the edge to the next undecimated site. With the help 
of (Al) it can be easily shown that an end site (— hoa^) 
decimation leads to the following flows: 

(f^y vKoD 

hpoi^if => (^) 2 hi Ji 2 h 2 a\al (A12) 

h J ih 1( T y a V i => \\J m ata z 2 . 
From the above equation we see that there is a factor 
(^r) associated with the displacement into the chain 



of the ctqUi , a-gcf operators. This is repeated partially 
in the case of a bond decimation of the — Joicrgcrf : 

hlJoi^o^l ^(01) J (01)20'foi) cr 2 

/ \ ( A1 3) 
h J 01 hi<jy<rl (^^J h {01) J (01)2 h 2 <T V {01) al. 



4. Leading Order Correlations 



Now that we know how the end operators neglected in 
the text evolve and discovered that their forms are invari- 
ant, we can show that these operators do not change the 
leading order contributions to the correlation functions 
of interest. The xx correlation function is 



C X L X =< h a x h L a x L > - < h a x >< h L a x L > 
= < G\(h aS- < h a x >){h L a x L - < h L a x L >)\G > 



= J2 <G\h a x \iP><iP\h L a x L \G>, 
with the sum over all excited states, ip. 



(A14) 



As the RG process progresses, all the end operators 
will be generated several times. However, we need only 
concern ourselves with the last set of these generated. 
Previously generated edge operators will have a larger 
suppression due to the more times they underwent edge 
displacement. 

Considering the last decimation step which is needed 
to obtain the x-x correlations and changing notation so 
that the remaining effective sites are I, r as in the text, 
we have: 
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**G •« •« (A15) 



The remaining low energy parts of the hamiltonian are 

H = -hiaf- Jitter* - h r a* . (A16) 

Two cases need to be considered: a site decimation and a 
bond decimation. In the case of a bond decimation, the 
ground state of the system is 

\H W > = -J=(| [t> | lr> +| ]l> | ]r>) (A17) 

and we immediately see that the only excitation that 
contributes in the sum (A15) is the afcr* term: 

C xx _ e -At-A rhgh ^ ( Alg ) 



which is the same result as was derived in the text while 
ignoring the additional edge operators from Eq. (A5). In 

h 2 

the above we used r^- <C 1. 

This demonstration can be repeated for the xx-z cor- 
relations and also carried to higher order with the same 
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